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Abstract 

The first LHC pp collisions at centre-of-mass energies of 0.9 and 2.36 TeV were 
recorded by the CMS detector in December 2009. The trajectories of charged particles 
produced in the collisions were reconstructed using the all-silicon Tracker and their 
momenta w^ere measured in the 3.8 T axial magnetic field. Results from the Tracker 
commissioning are presented including studies of timing, efficiency, signal-to-noise, 
resolution, and ionization energy. Reconstructed tracks are used to benchmark the 
performance in terms of track and vertex resolutions, reconstruction of decays, es- 
timation of ionization energy loss, as well as identification of photon conversions, 
nuclear interactions, and heavy-flavour decays. 
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1 Introduction 

The Compact Muon Solenoid (CMS) [l] is a general purpose detector at the Large Hadron 
Collider (LHC) of CERN. It has been designed primarily to perform new physics studies at 
the highest energies achievable with the LHC. The main components of CMS are a muon 
detection system, electromagnetic and hadronic calorimeters, and an inner tracking system 
(Tracker). The Tracker provides robust, efficient, and precise reconstruction of the charged 
particle trajectories inside a 3.8 T axial magnetic field. The nominal momentum resolution is 
typically 0.7 (5.0)% at 1 (1000) GeV/c in the central region and the impact parameter resolution 
for high-momentum tracks is typically 10 f/m. 

The reconstructed tracks of charged particles are among the most fundamental objects in the 
reconstruction of pp collisions. Tracks are used in the reconstruction of electrons, muons, 
hadrons, taus, and jets as well as in the determination of the primary interaction vertices. In 
addition, tracks may be used to identify b jets, in particular through evidence of a displaced 
vertex associated with a given jet. 

This paper describes the performance of the Tracker, which was evaluated with collision data 
from early LHC operations at centre-of-mass energies of 0.9 and 2.36 TeV. The next section 
contains a brief description of the Tracker. Section |3] illustrates the LHC data and conditions 
that underlie the analysis. Results obtained from the commissioning of the Pixel and Silicon 
Strip detectors are described in Section |4] Section |5] describes the track reconstruction and 
Section [6] presents tracking results demonstrating the overall performance of the Tracker. In 
particular, reconstructed tracks are used for track and vertex resolution measurements, the 
reconstruction of hadron decays, the estimation of ionization energy loss, the identification of 
photon conversions and nuclear interactions, and b tagging. Finally, conclusions are presented 
in Section [T] 

2 Tracker Description 

The CMS experiment uses a right-handed coordinate system, with the origin at the nominal 
interaction point, the x axis pointing to the centre of the LHC ring, the y axis pointing up 
(perpendicular to the LHC plane) and the z axis along the anticlockwise-beam direction. The 
azimuthal angle (p is measured in the xy plane, with cp = along the positive x axis and cp = 
n/2 along the positive y axis. 

The CMS Tracker [T], shown in Fig. Ill consists of two main detectors: a silicon pixel detector, 
covering the region from 4 to 15 cm in radius, and 49 cm on either side of the collision point 
along the LHC beam axis, and a silicon strip detector, covering the region from 25 to 110 cm in 
radius, and within 280 cm on either side of the collision point along the LHC beam axis. 

The CMS silicon pixel detector has 66 million active elements instrumenting a surface area of 
about 1 m^. It is designed to provide the determination of three high precision three-dimensional 
points on track trajectories. The detector consists of three concentric cylindrical barrel layers 
and four fan-blade disks which close the barrel ends. The barrel layers have an active length of 
53 cm and are located at average radii of 4.3, 7.3, and 10.2 cm. The endcap disks instrument the 
regions between radii 4.8 and 14.4 cm at mean longitudinal distances of 35.5 and 48.5 cm from 
the interaction point. The system provides efficient three-hit coverage in the region of pseudo- 
rapidity I //I < 2.2 and efficient tw^o-hit coverage in the region |f/| < 2.5. The active elements 
are n-in-n 100 //mx 150 jim pixels lIU w^hich are oriented with the smaller pitch in the azimuthal 
direction in the barrel and the radial direction in the disks. The 3.8 T magnetic field in CMS 
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Figure 1: r-z slice of the CMS Tracker. 



causes significant azimuthal Lorentz drift of the collected electrons in the pixel barrel which 
enhances the azimuthal charge sharing and therefore improves the resolution in that direction. 
The blades of the endcap disks are rotated by 20 degrees about their radial axes with respect to 
the disk planes to produce azimuthal charge sharing and radial Lorentz drift which enhances 
the radial charge sharing. The charge sharing improves the endcap resolution in both planes. 

The CMS silicon strip detector has 9.3 million active elements instrumenting a surface area of 
198 m-^. The detector consists of three large subsystems. The Tracker Inner Barrel and Disks 
(TIB/TID) extend in radius to 55 cm and are composed of four barrel layers, supplemented by 
three disks at each end. The TIB/TID delivers up to four r-(p measurements on a trajectory 
using 320 }im thick silicon microstrip sensors, which have their strips oriented parallel to the 
beam axis in the barrel and oriented radially in the disks. The strip pitch is 80 fim in the inner 
pair of TIB layers and 120 f/m in the outer pair of TIB layers. In the TID, the mean pitch varies 
between 100 f/m and 141 fim. The TIB/TID is enclosed within the Tracker Outer Barrel (TOB), 
w^hich has an outer radius of 116 cm. The TOB consists of six barrel layers of 500 f/m thick 
microstrip sensors w^ith strip pitches of 183 /^m in the first four layers and 122 jim in the last pair 
of layers. The TOB extends to ±118 cm in z. Beyond this z range, the Tracker EndCaps (TEC) 
instrument the region 124 < |z| < 280cm and 22.0 < r < 113.5cm. Each TEC is composed of 
nine disks that are instrumented with up to seven rings of radial-strip silicon detectors. The 
sensor thicknesses are thin (320 jim) in the inner four rings and thick (500 jim) in the outer three 
rings; the average radial strip pitch varies from 97 fim to 184 fim. The inner two layers of the 
TIB and TOB, the inner two rings of the TID and TEC, and the fifth ring of the TEC include a 
second microstrip detector module that is mounted back-to-back at a stereo angle of 100 mrad 
and that enables a measurement of the orthogonal coordinate. Assuming fully efficient planes 
and not counting hits in stereo modules, there are from 8 to 14 high precision measurements of 
track impact points for \rj\ < 2.4. 



3 Data Samples 

The results presented in this paper were obtained from data samples collected by the CMS 
experiment during LHC operation in December 2009 at proton-proton centre-of-mass energies 
of 0.9 and 2.36 TeV. The CMS axial magnetic field was maintained at the nominal value of 3.8 T 
and the silicon pixel and silicon strip detectors were biased at their nominal voltages. Due to 
the relatively low LHC luminosity, the CMS readout was triggered by the coincidence of signals 
from the beam scintillator counter (BSC) minimum bias trigger and the beam pick-up timing 
detector which detects the passage of the beam bunches [2J. The BSC minimum bias trigger 
requires that the arrival times of the signals from the forward and backward arms of the BSC 
(3.23 < \tj\ < 4.65) be consistent with the passage of particles emerging from a pp collision 
in the middle of CMS. In contrast, the BSC beam-gas trigger, used to veto non-collision events, 
requires that the arrival times be consistent with the passage of particles traversing the detector 
from one end to the other in time with particles from either beam. The total number of selected 
minimum bias events is approximately 305 000. 

Prior to the LHC pp collisions, the CMS experiment w^as commissioned using events containing 
cosmic muons during Cosmic Run At Four Tesla (CRAFT) [3J. The detector and magnetic field 
conditions during CRAFT were quite similar to the conditions during pp collisions. Thus, the 
results obtained from CRAFT provided good initial operating points for the pixel detector f4j, 
the strip detector [5], the tracker alignment [6J, and the magnetic field [7]. The data used in the 
referenced CRAFT papers were obtained in the fall of 2008, more than one year before the pp 
collisions. In most cases, more recent CRAFT data were used to improve on these results. 

4 Tracker Commissioning 

The follow^ing tw^o subsections describe the operating characteristics and performance of the 
silicon pixel and silicon strip detectors, respectively. 

4.1 Silicon Pixel Detector 

4.1.1 Operating Conditions 

In order to make maximal use of experience gained from the operation of the pixel detector 
•with cosmic rays during summer/autumn 2009, the operating conditions were not changed 
for the December 2009 data taking period. The coolant temperature w^as kept constant at 7°C. 
The bias potential applied to the 285 fim thick p-spray barrel sensors [8] was a uniform 150 V. 
The bias potential applied to the 270 }im thick p-stop endcap sensors [9J was a uniform 300 V. 
Small fractions of the barrel (1.0%) and endcap (3.1%) detectors were inactive resulting in a net 
operational fraction of 98.4% for the entire detector. 

The calibration procedures described in Ref. HJ were used to determine the ADC gains and 
pedestals for all channels. Iterative tuning reduced the mean (spread) of the readout threshold 
distributions for the pixel Readout Chips (ROCs) from the values measured during the 2008 
cosmic ray commissioning [4J to 2733 e (196 e) in the barrel detector and 2483 e (163 e) in the 
endcap detectors, w^here e is the magnitude of the electron charge. These measured threshold 
values apply only to the calibration procedure. Because the bandwidth of the preamplifiers is 
limited by pow^er considerations, small signals can take more than a bunch crossing time (25 
ns) to fire the zero-crossing discriminator that triggers the storage of the signal. This causes 
some smaller signals to be associated with the wrong bunch crossing and to be ignored by the 
readout system. The net result is that the effective or "in-time" thresholds are larger than the 
set values. 
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The effective thresholds are estimated by comparing the distribution of measured cluster x- 
sizes (azimuthal direction in the barrel detector and radial direction in the endcap detectors) 
with those predicted by the detailed pixel simulation, pixel AV lHOllTTl . The cluster sizes are 
sensitive to the effective thresholds. To avoid highly ionizing particles, the tracks used in this 
analysis were required to have momenta larger than 4 GeV/c. This selection ensures that even 
protons and deuterons produce signals that are w^ithin a few percent of the ionization min- 
imum. By varying the simulated thresholds until the measured and simulated distributions 
agree, the average effective thresholds are found to be approximately 3500 e in the barrel detec- 
tor and 3000 e in the endcap detectors. 

A study of the pixel hit reconstruction efficiency using a technique similar to the strip detector 



technique described in Section 4.2.4 suggests that the efficiency is larger than 99% for the live 
regions of the detector and is consistent with earlier work [2J. 



4.1 .2 Pixel Timing Scan 

The pixel detector readout system uses the 40 MHz LHC clock as input. Signals from the CMS 
trigger system must arrive at the correct time within the 25 ns clock cycle to associate the correct 
bunch crossing time stamp with any signal above the readout threshold. An optimally phased 
clock signal will maximize the number of pixels observed in clusters. The overall trigger timing 
was adjusted by varying the clock phase until the average barrel and endcap cluster sizes as 
measured in minimum bias triggers were maximized. These quantities are plotted versus clock 
phase in Fig. |2] The clock phase setting of 6 ns was found to optimize the smoothly varying 
detector averages. A finer module-by-module adjustment of the clock phase w^ill be performed 
when higher trigger rates become available. 
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Figure 2: The average cluster size distributions for the barrel and endcap pixel detectors in 
minimum bias events are plotted versus clock phase. 



4.1 .3 Operating Characteristics with Minimum Bias Triggers 

The distributions of the number of clusters observed in 0.9 TeV events selected by the minimum 
bias trigger are show^n in Figls] The observed data, shown as solid dots, are compared with fully 
simulated data, shown as histograms, that were generated with a recent tuning of the PYTHIA 
event generator [12]. The left plot shows the distribution for all events, whereas the right plot 
show^s the distribution after removing events that also satisfy the beam-gas trigger. There is an 
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excess of large multiplicity events that are removed by the beam-gas trigger requirement. The 
source of these events could be beam-gas interactions or beam scraping in the beam transport 
system near the interaction point. After removal of the beam background events, the measured 
distributions are approximately consistent with preliminary expectations. The measured aver- 
age cluster multiplicities per layer (barrel detector) and per disk (endcap detector) are listed in 
Table IT] They are compared with the expectation from the simulation and are found to be in 
rough agreement. It should be noted that the event generator is based on an event model that 
has not yet been tuned in detail and is not expected to provide accurate predictions. 
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Figure 3: The cluster multiplicity of (a) all minimum bias triggered events and (b) those that 
do not trigger the beam-gas veto in the 0.9 TeV data sample. The histograms show^ the similar 
distribution for a sample of simulated data. 



During the extremely low luminosity run in December 2009 (the instantaneous luminosity was 
typically in the range 10^^-10^'^ cm^-^s^^), the beam background events occurred at a rate that 
was roughly comparable to the rate of minimum bias triggers. Because they are characterized 
by particle trajectories that are nearly parallel to one of the beams, most background events 
(~90%) do not fire the minimum bias trigger but do have clusters in the endcap detectors 
and elongated clusters in the first tw^o layers of the barrel detector. At the beam energies of 
the December 2009 run, the pixel detector occupancies associated with the background events 
were typically five times larger than those associated w^ith minimum bias events. The beam- 
gas trigger veto effectively removes background events, as do cluster shape, track quality, and 
vertex requirements. 

The cluster charge distributions measured in the barrel and endcap detectors w^ith the 0.9 TeV 
sample are shown as solid dots in Fig. |4] Each entry is scaled by the ratio of the pixel sensor 
thickness to the track path length in the sensor. The solid histograms represent the expecta- 
tions from the PYTHlA-based, full detector simulation. The measured and simulated barrel 
distributions have similar peaks but the measured distribution is somew^hat broader than the 
simulated one. This may be due to residual pixel-to-pixel gain variation resulting from the use 
of a single gain for all 80 channels in each ROC column or residual module-to-module clock 
phase variation. The corresponding distributions for the endcap detectors have similar widths 
but indicate a 5% charge-scale mismatch. 
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Table 1: The average cluster multiplicity per layer/disk in 0.9 TeV minimum bias triggers. 
The simulation errors are entirely statistical and do not represent the uncertainties in the event 
modelling. The asymmetry seen in the forward and backward endcaps is caused by an offset 
in the luminous region along the beam axis. 
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Figure 4: The normalized cluster charge measured in the (a) barrel and (b) endcap pixel detec- 
tors for the sample of 0.9 TeV minimum bias events. The insets show the same distributions on 
semi-log scales. 
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4.1 .4 Lorentz Angle Calibration 

The use of n-in-n pixel technology and the large magnetic field in CMS imply that pixel hit 
reconstruction involves large Lorentz drift corrections (the typical bias corrections are 53 jim 
in the barrel and 10 /^m in the endcap). The estimation of track impact coordinates from pixel 
clusters is performed with two different algorithms. The simpler, faster "Generic Algorithm" 
||T3| uses the Lorentz width Wi to estimate the projected cluster size and bias correction. The 
Lorentz w^idth is the product of the effective thickness of the sensor Tgff and the tangent of the 
average Lorentz angle 6i: Wi = Tgff tan0£. Due to the focusing of the electric field at the n+ 
implants, the charge sharing near the n+ side of the sensors is reduced. This is modelled by the 
effective thickness which is 5-10% smaller than the physical thickness of the sensor substrate. 
The detailed PIXELAV simulation is used to extract the Lorentz w^idth by applying the Generic 
Algorithm to a sample of simulated clusters and by adjusting Wi to minimize the bias and 
maximize the resolution. The slower, more sophisticated "Template Algorithm" [14J fits pre- 
computed cluster shapes to the measured clusters. The Lorentz-drift effects are encoded in 
the cluster shapes and the same pixelav simulation is used to compute them. Therefore, the 
actual Lorentz calibration procedure is to tune the detailed simulation to agree with data and 
then to generate a Lorentz w^idth for the Generic Algorithm and cluster shapes for the Template 
Algorithm. 

Two different techniques have been used to perform the calibration. The 2008 cosmic ray data 



were calibrated by measuring the cluster x-sizes as functions of cot a (see Fig. 5a for definitions) 
and by determining the locations of the cluster-size minimum cotaniin[4]. In the pixel barrel, 
— cotamin is equal to tan^L = ^HfiB, where vh is the electron Hall factor, fi is the average 
electron mobility, and B is the magnetic field. The 2008 cosmic ray measurements suggested 
that the value of the electron Hall factor used in pixelav should be increased to 1.05 from the 
1.02 value determined in test beam measurements [15]. In 2009, the temperature of the detector 
was lowered and the bias voltage of the pixel barrel was increased, which changed the average 
Lorentz angles in both barrel and endcap detectors. New cosmic ray based determinations are 
reported in Table |2]and are compared with the tuned simulation. 

Table 2: The tangent of the Lorentz angle tan 9i as determined by 2009 calibrations. 

2009 Lorentz Angle Measurements 



Sample 


Detector 


Technique 


Measured tan 9i 


Simulation 


Cosmic Ray 


Barrel 


Cluster Size 


0.409±0.002(stat) 


0.407±0.002(stat) 


Cosmic Ray 


Endcap 


Cluster Size 


0.081 ±0.005(stat) 


0.080±0.004(stat) 


Minimum Bias 


Barrel 


Grazing Angle 


0.3985±0.0005(stat) 


0.4006±0.0005(stat) 


Minimum Bias 


Barrel 


Cluster Size 


0.409±0.002(stat) 


0.411±0.005(stat) 



The barrel calibration was repeated w^ith collision data in December 2009 using a new "grazing 
angle" technique [,16J . This technique makes use of the two-dimensional pixel segmentation to 
simultaneously measure the average transverse displacement of the charge carriers as a func- 
tion of distance along clusters produced by a sample of highly inclined tracks. Since longitu- 
dinal position in the cluster is completely correlated with depth in the junction, this technique 
determines the average transverse carrier displacement as a function of depth as shown graph- 
ically in Fig. |5b] The average Lorentz angle, extracted from the linear fit shown in the figure, 
is compared with the detailed simulation in Table [2| The extremely large population of highly 
curved, low transverse momentum tracks observed in minimum bias triggers spans the cot a. 
region needed to determine the minimum projected cluster size in the pixel barrel. This en- 
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ables the use of the cluster size technique as a cross check which is also reported in Table |2] 
Note that the two techniques are affected by different systematic effects and that a better than 
1% consistency is observed between the real and simulated measurements in all cases. A vari- 
ation of fitting procedures suggests that the total systematic uncertainty on the Lorentz angle 
calibration is less than 2%. 
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Figure 5: (a) The pixel local coordinate system and track angle definitions. The local z axis 
coincides with the sensor electric field E. The local x axis is chosen to be parallel to E x B where 
B is the axial magnetic field. The local y axis is defined to make a right-handed coordinate 
system. The angle a. is the angle between the x axis and the track projection on the local xz 
plane, (b) The transverse cluster displacement of highly inclined barrel clusters as a function 
of depth for a sample of 0.9 TeV minimum bias events at a magnetic field of 3.8 T The tangent 
of the Lorentz angle is given by the slope of a linear fit which is shown as the solid line. 
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4.1 .5 Resolution Study 

The intrinsic position resolution in a limited range of the angular acceptance was measured 
using tracks from minimum bias triggers that traverse overlapping sensors in the barrel layers. 
A similar analysis was performed in a very different angular region with 2008 cosmic ray data 
14] using the measurement technique given in Ref. I17| . Tracks passing through two overlap- 
ping modules in the same layer are used to compare the hit position with the expected position 
from the track trajectory. Because it is insensitive to alignment uncertainties, the difference of 
the local track impact points on a fitted trajectory is known about ten times more precisely 
than are the individual predicted hit positions. A double difference is formed by taking the 
difference between the measured hit position difference in the tw^o modules and the predicted 
trajectory position difference. The width of this double difference distribution is insensitive to 
translational misalignment of the overlapping modules. 

To limit the effect of multiple scattering, a minimum track momentum of 2.5 GeV/c is required. 
Clusters with measured charge below 10 000 e or containing pixels on the sensor edges are 
excluded. The double difference widths are fitted with a Gaussian and the uncertainty from 
the trajectory prediction is subtracted quadratically to recover the hit resolution on the position 
difference. With the assumption of equal resolution for each of the modules in the overlap, 
the final fit values for the resolution for a single module are 12.8 ± 0.9 }im along x and 32.4 ± 
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1.4 f/m along y. The pixelav simulation is used to generate a sample of clusters that has the 
same distribution of impact angles as the measured sample. Since the simulation does not 
include the double-size pixels that span the gaps between the sixteen readout chips which tile 
each module, a subsample of the overlap data sample is used to determine single-size-pixel 
resolutions of 12.7 ± 2.3 f/m along x and 28.2 ±1.9 ^m along y. These numbers can be directly 
compared with those extracted from Gaussian fits to the simulated residual distributions. The 
simulated resolutions are 14.1 ± 0.5 f/m and 24.1 ± 0.5 jim along x and y, respectively, and agree 
reasonably well with the measured resolutions. Because overlaps occur only at the edges of the 
track a-angle acceptance w^here the x sizes of the clusters deviate from the optimal size of tw^o, 
the measured and simulated x resolutions are somewhat w^orse than the typical x resolution 
(less than 10 jim) expected for most collision-related clusters. The measured and simulated y 
resolutions are expected to be typical of the detector performance. 

4.2 Silicon Strip Detector 

4.2.1 Operating Conditions 

All of the modules in the strip tracker were biased at 300 V in the early collision running. This is 
the same setting that w^as used in the CRAFT studies and is well above the full depletion voltage 
for the sensors. Similarly, the coolant temperature was set at 4-6 °C, the same as in the CRAFT 
study. This meant that the p+ on n sensors IITSi were approximately at room temperature. 

As described in the Technical Design Report for the CMS Tracker Ill9l |20| , there are two main 
modes of operation for the strip tracker analogue pipeline integrated circuits (APV25 ||2T|): 
peak and deconvolution. In deconvolution mode, the output charge for each strip represents 
a weighted sum of three consecutive pipeline cells [22 J. Although deconvolution mode was 
designed to avoid signal pile-up in high (design) luminosity operations, it will be necessary to 
run in this mode whenever the expected separation betw^een proton bunches will be less than 
a few^ hundred nanoseconds. The luminosity in the early collision running was very low^ and 
the bunches well separated; most of the strip data were collected in peak mode, which is based 
on the signal in a single pipeline cell. All of the data, whether in peak or deconvolution mode, 
were zero suppressed, meaning that only strips w^hich were part of clusters were read out for 
each event. 

Many of the starting parameters for the strip tracker during the early collision running had 
been established in the preceding CRAFT period. For example, the timing of the tracker sub- 
systems (in peak mode) w^ith respect to CMS triggers w^as set during the cosmic-ray muon 
studies. Similarly, the alignment parameters for the strip detector modules were derived from 
the same studies. 

As part of the alignment process, offsets had been determined for cluster positions in sensors 
due to the Lorentz drift of holes and electrons under the influence of the solenoid field. For 
barrel layers, the Lorentz angle correction for cluster positions during track reconstruction is 
about 10 f/m, which is significantly larger than the 3-4 jim alignment precision achieved in the 
cosmic ray studies [6J. 

4.2.2 Strip Timing Scan 

As the strip tracker was operated in peak mode at the start of the early collision running, the 
trigger timing established in the preceding CRAFT period could be used. In CRAFT the sam- 
pling time of the APV25's was set within each subsystem by means of a dedicated synchro- 
nization signal, adjusted according to the measured readout fibre lengths. The synchronization 
of the subsystems was obtained using the signal from cosmic-ray muon tracks. Details on how 
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the scan was done can be found in Ref. Il23l . Toward the end of the data collection period the 
APV25 mode was changed from peak to deconvolution and since timing is more critical in the 
latter, a fine-delay scan was made following the mode change. For expediency only one layer 
(TOB L3) was used in the study. 

Figure [6] shows the result of the fine-delay timing scan. The timing adjustment for the clock 
and trigger signals is set on the front-end hybrids and the smallest step size is 1.04 ns. From 
the figure it can be seen that the timing prior to the scan had been off by about 10 ns from 
ideal. This level of mistiming resulted in an estimated 2.5% decrease in Signal-to-Noise (S/N) 
in the strip modules during the peak mode running, where the delay timing is less critical. The 
amplitude that is measured in the timing scan represents the signal of the highest pulse height 
strip in a cluster scaled by the ratio of the sensor thickness to the path length of the track in the 
sensor. 

Following the scan, the timing offsets for all of the strip tracker subsystems w^ere updated and 
some data were collected in deconvolution mode. No data samples were collected in peak 
mode with the new delays. 
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Figure 6: Normal-incidence-scaled charge (arbitrary units) of the highest pulse height strip in 
a cluster as a function of the readout delay with respect to the CMS trigger, in deconvolution 
mode. The dashed vertical line corresponds to the setting prior to the timing scan. 



4.2.3 Signal-to-Noise Measurements 

Signal-to-Noise measurements w^ere made in both peak and deconvolution modes. In peak 
mode, the S/N ratio was determined at both centre-of-mass energies, 0.9 and 2.36 TeV, whereas 
deconvolution mode is restricted to 2.36 TeV. The ratio is evaluated on the basis of charge clus- 
ters associated with reconstructed tracks, where the individual strip noise values are taken 
from calibration runs. For track angles that are not normal to the surface of modules the signal 
values are scaled by the cosine of the angle relative to the local normal. This is done to give the 
same expectation value per cluster for modules of the same type. Cluster noise, which takes 
into account the noise of each strip w^ithin a cluster, is used as the denominator in the S/N ratio. 
When all strips within a cluster have the same noise, cluster noise is equivalent to the noise of 
a single strip. Further details on the determination of the S/N ratio can be found in Ref. ||24|. 



Figures 7a and 7b show the S/N distributions for the TIB and TOB modules, respectively, in 
deconvolution mode. Included with each distribution is the result of the fit to a Landau distri- 
bution convolved with a Gaussian distribution. The most probable value of the fitted curves is 
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taken to be the S/N value and results for all of the strip tracker subdetectors are summarized in 
Table [Slfor all three running conditions. Peak values shown in the table have not been corrected 
for the 2.5% loss due to non-optimal timing. They are comparable with results obtained in the 
CRAFT studies and in earlier cosmic ray studies. The difference in peak and deconvolution 
mode S/N values stems largely from the higher noise in deconvolution. After calibration there 
is some variation in signal values (measured in electrons) for the two modes, but this has been 
shown to be within 10%. The S/N ratio should not depend on the centre-of-mass energy and 
this is confirmed by the table entries. 

Although it is not possible to directly compare channel noise distributions in the early collision 
data with results from calibration runs given the zero suppression, the frequency and distribu- 
tion of clusters in empty LHC buckets provide an indirect cross-check of the calibration results 
and assumptions about the Gaussian and uncorrelated nature of the noise. For example, with 
bad modules excluded from the readout the mean number of clusters in empty buckets, out 
of some 9 million channels, w^as 4.2. This is consistent w^ith the clustering rules, w^hich require 
a certain number of standard deviations (five for the total charge in a cluster), and Gaussian 
probabilities. By way of contrast, there were ~1200 clusters per minimum bias trigger in the 
0.9 TeV data. 
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Figure 7: Signal-to-Noise distributions in deconvolution mode for (a) (thin sensor) TIB and 
(b) (thick sensor) TOB modules. The curves are results of the fits to a Landau distribution 
convoluted with a Gaussian distribution. 



Table 3: Summary of strip tracker Signal-to-Noise measurements. The peak mode ratios have 
not been corrected for the estimated 2.5% decrease in signal from the trigger mistiming, as 
described in the text. 



Conditions 



TIB TID TOB TEC thin TEC thick 



0.9 TeV, peak mode 27.4 26.7 34.1 28.8 35.7 

2.36 TeV, peak mode 27.4 26.8 34.1 28.8 35.7 

2.36 TeV, deco mode 20.3 19.2 23.9 20.3 26.1 



4.2.4 Strip Layer Efficiencies 

Efficiencies for strip tracker layers were determined using events that were collected in peak 
mode. Reconstructed tracks in these events were required to have a minimum of 8 hits in order 
to be used in the efficiency measurements. To avoid inactive regions and allow for alignment 
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imprecision, trajectories passing near the edges of sensors were excluded. The presence of a 
hit anywhere within the non-excluded region of a traversed module was counted as a posi- 
tive response; efficiency is determined by the ratio of positive responses to the total number 
of traversing tracks. Layers under study were not removed from the track reconstruction and 
could in fact count toward the minimum hit requirement. The total integrated hit efficiency 
during the early collision period w^as measured to be 97.8%, which is essentially explained by 
the number of bad modules in the strip tracker. That is, about 2.2% of the modules have been 
excluded from the readout because of problems with high voltage short circuits, control ring 
failures, or other issues. With known problem modules excluded, the overall hit efficiency is 
99.8%, consistent with the ~0.2% bad channel rate from the construction process. Detailed sim- 
ulations, used to determine track reconstruction efficiency, take into account inactive regions 
in addition to the measured efficiencies. The efficiency measurements for the collision data 
include an estimated 0.04% systematic error due to the use of the layers under study in the 
reconstruction process and the wide search w^indows w^ithin modules. 

4.2.5 Energy Loss Measurement 

Although the primary function of the strip tracker is to provide hit position information for 
track reconstruction and precise momentum determination, the wide linear range of the strip 
channel output also provides a measure of energy loss. That is, the charge collected in a hit 
cluster is directly proportional to energy lost by a particle, largely through ionization, while 
traversing the silicon. For reconstructed tracks the angle 9 between the track direction and 
the axis normal to module sensor is w^ell defined for each hit on the track. The instantaneous 
energy loss per unit path length (dE/dx) in the silicon is then approximated by the quantity 
A£/ (AL • sec 6), where A£ is the cluster charge expressed in units of MeV and AL is the normal- 
angle thickness of the active volume of the silicon sensor. All of the TIB and TID modules 
and the modules on rings 1-4 of the TEC have silicon sensors that are 320 f/m thick, whereas 
the TOB and TEC ring 5-7 modules have 500 }im thick sensors. Some 30 fim of the nominal 
thicknesses for both thin and thick types is inactive material, i.e., does not contribute to the 
charge collection. 

In zero-suppressed readout, which was used exclusively in the early collision period, there 
are 8 ADC bits for the charge on each channel within a cluster. Channel gains are set such 
that a single ADC count corresponds to about one-quarter of the average noise and full scale 
corresponds to approximately three times the average loss expected from normally incident 
minimum ionizing particles. The highest two ADC values have a special significance: 254 
implies a value between 254 and 1024 counts, and 255 indicates that the actual value was in 
excess of 1024 counts. The dE/dx algorithm includes the saturated values but w^ithout any 
special treatment. 

The main point in determining energy loss per unit path length is that, for a given medium, 
dE/dx depends largely on the velocity (j6) of the traversing particle. By combining dE/dx 
information with the measured momentum p of a track, one can determine the mass of the 
traversing particle. On the scale of charged particle momenta in CMS collisions, there is only 
a limited range near the low end where the difference in /S values is significant enough to dis- 
tinguish among long-lived hadrons. The momentum range w^here pions would have relatively 
large energy loss is such that tracks tend to curl up in the 3.8 T solenoid field and thus fail to be 
reconstructed. 

The strip hits on reconstructed tracks represent independent measures of dE/dx, ignoring the 
negligible loss of energy in traversing the tracker. Although pixel hits are included in the track 
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reconstruction, they are not used in the dE/dx calculation due to their more limited linear 
range. Several methods have been used to determine an estimate for the most probable dE/dx 



value based on the measurements in the strip tracker modules traversed by a track. Figure 8a 
for example, shows the relationship between the Harmonic-2 dE/dx estimator |25| and mo- 
mentum for 0.9 TeV data taken in peak mode. In the figure, clear bands can be seen for kaons 
and protons and to a much lesser extent for deuterons. 

An estimate of the mass of each candidate can be obtained using the particle momentum and 
the measurement of the ionization energy loss provided by the dE/dx estimators. To this 
end the following relation betw^een dE/dx, p, and m is assumed for the momenta below^ the 
minimum-ionizing region: 



dE 
dx 



K- 



m^ 



C 



(1) 



The proton line in Fig. 8a is used to extract the parameters K and C in Eq.lT] The 0.7-1.0 GeV/c 
range in the proton band is used for the reference data fit, while extrapolations based on the 
same K and C values yield a good agreement for protons with momenta above and below the 
reference range and for kaons. 

The mass spectrum that results from inverting Eq.lTlfor all tracks with dE/dx > 4.15MeV/cm 
and p < 2 GeV/c is shown in Fig. 8b From the frequency plot one can observe clear kaon and 
proton peaks as well as good agreement for the peaks from a Monte Carlo simulation. There is 
also evidence for a deuteron peak in data, although saturation of the ADC scale is particularly 
pronounced for deuterons given their reduced /3 values and relatively higher | rj \ values. That 
the deuteron peak is poorly modelled by the simulation is partly understood as the underlying 
generator, PYTHIA, does not produce deuterons by design, although they can be produced in 
the subsequent GEANT [26J hadron show^ers. 
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Figure 8: Energy loss versus the momentum of tracks (a) and frequency of tracks as a function 
of track mass as determined from the measured energy loss and momentum (b). The lightly 
shaded line in the (a) indicates the fit in the reference range of the proton band while the darker 
lines correspond to extrapolations for kaons, protons, and deuterons based on the fit parame- 
ters. 
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5 Track Reconstruction 

The track reconstruction algorithms rely on a good estimate of the proton-proton interaction 
region, referred to as the beamspot. The beamspot is used as a precise estimate of the primary 
interaction point (in the transverse direction) prior to primary vertex reconstruction and as 
the sole primary interaction point if no primary vertex is found. When the beamspot centre 
is displaced from the expected position there is a correlation between the transverse impact 
parameter (dxy) and the angle of the track at the point of closest approach {(po). The beamspot 
fitter [27j uses an iterative x^ fitter to exploit this correlation between dxy and (po, looping over 
a sample of reconstructed tracks (using the old beamspot) to determine the new beamspot 
parameters. After the beamspot is measured, the standard track reconstruction is performed. 
During the 2009 data-taking, a beamspot was fitted during each LHC fill; fill-to-fill variations 
were at the level of ~0.5 mm in x and y, and ~2 cm in z. 

Starting from the location of the beamspot, an initial round of track and vertex reconstruction is 
performed using only pixel hits. The pixel vertices found at this stage are used in the standard 
tracking. The standard track reconstruction at CMS is performed by the combinatorial track 
finder (CTF) Il28l . Tracks are seeded from either triplets of hits in the tracker or pairs of hits 
with an additional constraint from the beamspot or a pixel vertex, yielding an initial estimate 
of the trajectory, including its uncertainty. The seed is then propagated outward in a search for 
compatible hits. As hits are found, they are added to the trajectory and the track parameters 
and uncertainties are updated. This search continues until either the boundary of the tracker 
is reached or no more compatible hits can be found. An additional search for hits is performed 
starting from the outermost hits and propagating inward. In the final step, the collection of hits 
is fit to obtain the best estimate of the track parameters. 

The current implementation of the CTF performs six iterations. Between each iteration, hits 
that can be unambiguously assigned to tracks in the previous iteration are removed from the 
collection of tracker hits to create a smaller collection that can be used in the subsequent iter- 
ation. At the end of each iteration, the reconstructed tracks are filtered to remove tracks that 
are likely fakes and to provide a means of quantifying the quality of the remaining tracks. The 
filtering uses information on the number of hits, the normalized x^ of the track, and the com- 
patibility of the track originating from a pixel vertex. Tracks that pass the tightest selection are 
labelled highPurity. The first two iterations use pixel triplets and pixel pairs as seeds to find 
prompt tracks with pj > 0.9 GeV/c. The next iteration uses pixel triplet seeds to reconstruct 
low-momentum prompt tracks. The follow^ing iteration uses combinations of pixel and strip 
layers as seeds, and is primarily intended to find displaced tracks. The final two iterations use 
seeds of strip pairs to reconstruct tracks lacking pixel hits. 

6 Tracking Performance 

The results presented here come from the sample described in Section |3| using data taken at 
both centre-of-mass energies (0.9 and 2.36 TeV), unless stated otherwise. To reduce the back- 



ground from beam-gas events, discussed in Section 4.1.3 and to select useful events for track- 
ing studies, two additional criteria are imposed for most of the results in this section. First, 
more than 20% of the reconstructed tracks in an event must be flagged as highPurity if there are 
at least 10 tracks in the event. Second, a primary vertex must be reconstructed in the region of 
pp interactions (see Section [0]|. 



The alignment parameters for the Tracker were computed from approximately two million 
cosmic ray muon tracks collected during CRAFT running in November 2009 as described in 
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Section [3] The nominal values of the alignment parameter errors have been used in the track 
reconstruction. Since the applied procedure was similar to the one discussed in Ref. HI, the 
resulting precision is also very similar. In particular, the width of the distribution of the mean 
of the residuals (taken as a measure of the local alignment precision) in the pixel barrel local x 
and y coordinates is 3 /^m and 4 ^m, respectively. 

The simulated events are minimum-bias events produced w^ith the pythia 6.4 1129 1 event gen- 
erator, tune D6T 1130 1 , at centre-of-mass energies of 0.9 and 2.36 TeV (10 million events each) 
and processed w^ith a simulation of the CMS detector response based on GEANT4. The mis- 
alignment, miscalibration, and dead-channel map corresponding to the detector status and 
calibration accuracy at the time of the first LHC collisions have been included in the simula- 
tion. The longitudinal distribution of the primary collision vertices has been adjusted to match 
the data. 

6.1 Basic Tracking Distributions 

The highPurity tracks are selected, w^ith additional requirements of |(iz| < 10 cr^ (where dz is 
the longitudinal impact parameter with respect to the primary vertex and Hz is the combined 
track and primary vertex uncertainty in z) and (ipj / px < 10%, to compare the data and sim- 
ulation. Figure [9] shows the results of this comparison for several important track parameters. 



The distribution of the number of tracks per event, shown in Fig. 9a has been normalized to 
the number of events. The data clearly have more tracks per event than are present in the sim- 
ulated data. This is believed to be due to an as-yet unoptimized tune of the pythia generator. 
To be able to compare shapes, the other distributions have been normalized to the number of 
reconstructed tracks in the data. There is general agreement betw^een the data and simulation 
distribution shapes for all other tracking variables. In particular, the features in the (p distribu- 
tion, due to inactive modules, are well modelled by the simulation. 

6.2 Primary Vertex Resolution 

The reconstruction of the primary interaction vertex in the event starts from the track collection. 
The tracks are clustered based on the z coordinate of the track at the point of closest approach 
to the beamline. The clusters are fit with an adaptive vertex fit [31J, where tracks in the vertex 
are assigned a weight between and 1 based on their proximity to the common vertex. 

The primary vertex resolution strongly depends on the number of tracks used in fitting the 
vertex and on their pj. To measure the resolution, the tracks in an event w^ith only one vertex 
are randomly split into two different sets and used to independently fit the primary vertex. 
The distribution of the difference in the fitted vertex positions can then be used to extract the 
resolution by fitting a Gaussian to it and dividing a by v 2. To examine the effect of the pj 
of the tracks in the vertex, we study the resolution versus the number of tracks in the vertex 



for different average pj of tracks in the vertex. Figure 10 shows the x, y, and z resolutions for 
different average pj ranges. While the resolution differs considerably depending on pj and 
multiplicity, the simulation accurately reproduces the data results. 

6.3 Reconstruction of Particle Decays 

6.3.1 V Reconstruction 

V*^ particles are long-lived (ct > 1 cm) neutral particles reconstructed by their decay to two 
charged particleaM Kg — ?► Tz^n^ and A*^ — )■ prc^ . Reconstruction of V*^ decays requires finding 

^Charge conjugate states are implied throughout the paper. 
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Figure 9: Comparison of the data (points) and simulation (histogram) distributions of tracking 
parameters: (a) number of tracks per event, (b) number of hits used per track, transverse (c) 
momentum pj, (d) track pseudorapidity rj, (e) azimuthal angle (p, (f) transverse impact param- 
eter d-xy with respect to the primary vertex, (g) longitudinal impact parameter d^ with respect 
to the primary vertex, and (h) normalized x^. The simulated distributions are normalized by 
area to the data distributions. 
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Figure 10: Primary vertex resolution distributions in (a) x, (b) y, and (c) z versus number of 
tracks. The three sets of results in each plot show different average pj ranges and within each 
Pt range, data and simulation are compared. 



oppositely charged tracks that are detached from the primary vertex and form a good sec- 
ondary vertex w^ith an appropriate invariant mass. For the A^, the lowest momentum track is 
assumed to be the pion. As no further particle identification is required, a V° candidate can 
appear in both Kg and A° samples. To be considered as a V*^ decay track, a track must have at 
least 6 hits, a normalized x^ less than 5, and a transverse impact parameter with respect to the 
beamspot greater than 0.5(Tip, where ajp is the calculated uncertainty (including beamspot and 
track uncertainties). The reconstructed V^ decay vertex must have a normalized x^ less than 
7 and a transverse separation from the beamspot greater than 15(Tt, where crj is the calculated 
uncertainty (including beamspot and vertex uncertainties). In addition, the V'^ candidate is 
discarded if either of the daughter tracks has hits that are more than 4(73d from the V'^ vertex, 
towards the primary vertex, where ct^d is the uncertainty in the vertex position. 

The mass resolution of the V*^ depends on t] as well as on the decay vertex position and a sin- 
gle Gaussian is not a sufficiently accurate functional form for the signal. Therefore, a double 
Gaussian with the same mean was used to fit the signal. For the background shapes, a linear 



background was used for 7r+ n and the function a I 



m — m, 



nij 



' was used for the pre spec- 



trum where m is the pn^ invariant mass and a and h are free parameters. The W^n^ and pn^ 
mass distributions, along with the overlaid fits, are shown in Figs. 11a and llb[ respectively. 
Tables m and m show the reconstructed V*^ masses and resolutions obtained from the data and 
simulation. While the various results are close to expectations, significant discrepancies are 
present. These features can be examined as a function of track kinematic variables to better 
understand the CMS tracker and magnetic field. This work is ongoing. 

Table 4: Masses obtained from data, world average [32 1, and simulation (reconstructed and 
generated). The uncertainties for data and simulation results are statistical only. 



VO 


Data 


Mass (Me^ 
PDG 


Simulation 


Generated 


K^ 
AO 


497.68 ± 0.06 
1115.97 ±0.06 


497.61 ± 0.02 
1115.683 ±0.006 


498.11 ±0.01 
1115.93 ±0.02 


497.670 
1115.680 
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Figure 11: The invariant mass distributions of (a) tz^tz with a fit to the Kg and (b) pn with a 
fit to the A°. 



Table 5: V^ mass resolutions obtained from data and simulation. The narrow and wide Gaus- 
sian resolutions are Ci and (72, respectively. The di fraction is the fraction of the yield from the 
narrow Gaussian. The final row gives the average resolution, obtained from the square root of 
the weighted average of the two resolutions squared. Uncertainties are statistical only. 



Parameter 


K^ Data 


Kg Simulation 


A" Data 


A*^ Simulation 


ai{MeV/c^) 


4.53 ±0.12 


4.47 ± 0.04 


1.00 ±0.26 


1.71 ± 0.05 


t7-2(MeV/c2) 


11.09 ±0.41 


10.49 ± 0.11 


3.25 ±0.14 


3.71 ± 0.09 


(Ti fraction 


0.58 ± 0.03 


0.58 ± 0.01 


0.15 ±0.05 


0.44 ± 0.03 


^(MeV/c^) 


7.99 ± 0.14 


7.63 ± 0.03 


3.01 ± 0.08 


2.99 ± 0.03 



6.3.2 V Lifetime 

For the 0.9 TeV centre-of-mass energy data and simulation, invariant mass distributions are 
made for different bins of proper decay length, ct = mcL/p, where L is the measured decay 
length. These distributions are fitted to obtain the yield, leading to the uncorrected ct distribu- 
tion as seen in Fig. 12a for the Kg data. The uncorrected ct distribution from the simulation is 
divided by the generated exponential shape given by e^'^^/'^^s/m to obtain the correction factor 
versus ct. The uncorrected data ct distribution is divided by the correction factor to obtain the 

for the Kg. This distribution is fitted with an expo- 



12b 



corrected ct distribution as seen in Fig. 
nential, the slope of which gives the measured lifetime. The good fit to an exponential func- 
tion (;i(;^/NDOF = 8.1/8) indicates that the simulation accurately reproduces the efficiency 
variation versus lifetime. The fitted results, Tj^o = 90.0 ± 2.1 ps and t^o = 271 ± 20 ps (with 

X^/NDOV = 11.3/6), are both within 1 C7 of the world average ||32l. 



6.3.3 Reconstruction of K*(892)- and E" 

The reconstructed sample of V° particles was exploited to reconstruct decays of other particles. 

The Kg candidates are combined with charged tracks from the primary vertex to search for the 



strong decay K* (892) 



Kin-. 



For this analysis, events were required to contain a recon- 
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Figure 12: Kg ct distributions for (a) uncorrected data and (b) corrected data with an exponen- 
tial fit. 



structed primary vertex consisting of more than two tracks and a fit pr obabi lity greater than 
0.5%. The Kg candidate must pass the same criteria as described in Sec. 6.3.1 In addition, the 



requirement on the impact parameter significance of the pions from the Kg is increased from 0.5 
to 2. The Kg candidates must also have a mass within 20 MeV/c'^ of the nominal mass and the 
Kg flight path must pass within 2 mm of the primary vertex. The charged track in the K* (892)^ 
decay must have a normalized x^ less than 2, at least two hits in the pixel detector, at least 
seven total hits, pj > 0-5 GeV/c, \rj\ < 2, and pass within 2 (3) mm of the primary vertex in the 
direction transverse to (along) the beam line. The Kg/r^ invari ant m ass is calculated using the 

The figure also shows an 
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world-average value of the Kg mass |32| and is shown in Fig. 

overlay of a fit to the Kg tt^ mass distribution. The fit uses a Breit-Wigner for the signal plus a 

threshold function for the background 



m^ 



Ml 



+n 



kM* 



+ B 



exp 



MK + Mn-m 



where m is the Kg/r^ invariant mass, Mk* and F^* are the mass and width of the K*(892)^, 
Mk and M^ are the world-average masses of K^ and n^ , and S, B, and p are free parameters. 
The K* width (Fj<:* ) is fixed at the w^orld average value of 50.8 MeV/c^ ||32| , w^hile the K* mass 
{Mk* ) is a free parameter. The mass returned by the fit, 888.3 ± 3.2 MeV/c^, is consistent with 
the world average value of 891.66 ± 0.26 MeV/c^ 1521. 

The E^ was reconstructed through its decay to A^tt^. The E^ is a long-lived baryon, with 
a decay topology different from that of the K*(892)^: the n^ from the E^ decay should be 
detached from the primary vertex rather than originating from it. The A*^ candidates were 



reconstructed as described in Sec. 6.3.1 except that a looser transverse significance cut of 10 
(rather than 15) w^as applied. A*^ candidates with a mass within 8 MeV/c^ of the world-average 
value were combined with charged tracks with the same sign as the pion in the A*^ decay. The 
A^n^ fit used a A^ mass constraint and the vertex was required to have a fit probability better 
than 1%. All three tracks involved in the decay w^ere required to have at least 6 valid hits and 
a 3D impact parameter with respect to the primary vertex greater than 3cr. The resulting mass 



plot, shown in Fig. 13b is fit with a single Gaussian for the signal and a background shape of 
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Aq^^^^^ + Bq^^^'^'> where q = m — Ma — M^, m is the A^/r^ invariant mass, and A and B are free 
parameters. The measured mass of 1322.8 ± 0.8 MeV/c^ is close to the world average value of 
1321.71 ± O.OZMeV/c^ [32J. The resolution of 4.0 ± 0.8MeV/c2 is consistent with the simulation 
result of 3.6 ± 0.4 MeV/c2. 
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Figure 13: Invariant mass plots of (a) Kg/r with a fit to the K*(892) and (b) A^k with a fit 
to the H". 



6.4 Particle Identification Using Measured Energy Losses 

Estimating the energy loss (dE/dx) of a particle by means of charge collected by the CMS sili- 
con strip tracker is described in Sec. [4.2.5 In this section, applications of dE/dx measurements 
are used to identify protons and kaons produced in A*^ and (p decays. 



6.4.1 dE/dx Verification with A^- prt Decays 

The kinematics of the A^ — > pn^ decay requires pp > p^ for all A*^ particles reconstructed 
at CMS. This provides a clean source of protons and pions which can be used to check the 



dE/dx results. We apply the same selection as in Section 6.3.1 and plot the dE/dx distribution 
as a function of the momentum for tracks associated to V^ ca ndid ates in the mass range 1.08- 
1.16 GeV/c^, separately for the highest momentum tracks (Fig. 14a i and the lowest momentum 
tracks (Fig. 14b| . As expected, the highest momentum tracks are generally found near the 
proton curve while the low^est momentum tracks are generally inconsistent with the proton 
curve. The few^ exceptions are consistent with background under the A*^ peak. 

6.4.2 Reconstruction of ^(1020) -^ K+K" 

The ^(1020) — > K+K^ decay w^as reconstructed using data taken at 0.9 TeV centre-of-mass en- 
ergy. The candidate kaon tracks come from the collection of highPurity tracks and are required 
to have pj > 0.5GeV/c, normalized X'^ < 2, at least five hits, \ri\ < 2, and a transverse impact 
parameter with respect to the reconstructed beamspot smaller than 3 mm. Finally, for tracks 
•with p < 1 GeV/c, the track must have a measured dE/dx consistent with the kaon hypothesis 
(see Eq.|l): KiM'^^/p)^ + C < dE/dx < K{M^''''/pY + C. The parameters of the dE/dx cut 



for kaons are those extracted from a fit to the dE/dx vs. p distribution, as described in Sec. 4.2.5 



6.5 Reconstruction of Photon Conversions and Nuclear Interactions 
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Figure 14: Estimated energy loss as a function of the momentum for (a) the highest momentum 
track and (b) the lowest momentum track for the A*^ candidate decay products. The superim- 



posed curves comes from the proton fit in the inclusive track sample shown in Fig. 8a 



We use a compatibility window of ±200 MeV/c^ around the K mass, with M^"^ and M'"^" being 
lower and upper boundaries of this window. 

The fit of the mass spectra of pairs of tracks accepted by the dE/dx selection used the sum 
of two normalized functions: a convolution of a relativistic Breit-Wigner shape with a Gaus- 
sian for the (p signal and a n arc tangent function for the background. The mass plot and over- 
laid fit are shown in Fig. 15a The fitted (p mass of 1019.58 ± 0.22 MeV/c-^ is in agreement 
•with, the w^orld-average value of 1019.455 ± 0.020 MeV/c^. The resolution found in data is 
1.29 ± 0.32 MeV/c^, in agreement with the value found in simulation, 1.41 MeV/c^. Candi- 



dates in which at least one track fails the dE/dx requirement are shown in Fig. 15b where 
only background is observed, indicating that the dE/dx requirement has a high efficiency to 
select (/>(1020) candidates. 

6.5 Reconstruction of Photon Conversions and Nuclear Interactions 

While the tracker is essential for finding charged particles and measuring their momenta, the 
tracker material is also a source for interactions. For photons, interactions with the tracker 
material can produce e+e^ conversion pairs, while for hadrons, nuclear interactions can pro- 
duce multiple hadrons. Photon conversions in the Tracker reduce the efficiency for low-energy- 
photon finding by the electromagnetic calorimeter, while nuclear interactions reduce track find- 
ing efficiency and can affect the resolution of many hadronic observables such as jets or missing 
transverse energy. Thus, identification of conversions and nuclear interactions can be used to 
improve many aspects of the event reconstruction. Furthermore, studies of conversions and 
interactions can be used to improve our understanding of the material in the Tracker. 

The electrons and positrons from converted photons can be identified by the electromagnetic 
calorimeter and used as seeds for track reconstruction 1.33.1. In the minimum bias events col- 



lected in December 2009, however, the photons have a soft spectrum as seen in Fig. 16a and 
therefore the conversion pairs are unlikely to reach the electromagnetic calorimeter. These 
conversion pairs can still be reconstructed by using tracker-seeded conversion reconstruction 
techniques, made possible by the iterative tracking algorithm described in Section [5] which 
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Figure 15: K+K^ invariant mass distribution, with (a) both kaons satisfying the dE/dx require- 
ment and with (b) at least one particle failing that requirement. In (a) a fit to the ^(1020) 
hypothesis is show^n. 



extends the capability of reconstructing low-px and detached tracks. The essential signature 
of a massless conversion photon is the two parallel tracks at the production vertex, in both 



the transverse and longitudinal planes. The reconstructed invariant mass, shown in Fig. 16b 



shows the effect of the mass resolution, which is well modelled by the simulation. Two dif- 
ferent conversion reconstruction approaches have been used. Both methods fit two oppositely 
charged tracks to a common 3D vertex with the constraint that the two tracks are parallel at the 
vertex. The met hod s diff er mainly in the preselection of the track pairs. The first method, from 
which Figs. 16a and 16b are derived, requires both tracks have at least 3 hits and normalized x^ 



less than 10 and at least one track with 5 or more hits. The tracks are required to have positive 
charge-signed transverse impact parameter, positive distance of minimum approach in 2D (i.e., 
the two full track circles have one or no intersection in the transverse plane), small z separation 
at their innermost point (|Az| < 5 cm) if they are in the barrel, and a small opening angle in 
both the transverse (Ai^ < 0.2) and longitudinal plane (Acot0 < 0.1 where 6 is the polar angle 
relative to the z axis). The vertex fit must have a x^ probability better than 5 x 10^^ and be 
located inside the innermost hits on the tracks. To increase efficiency, the second method takes 
all tracks with a X'^ probability above 10^^ and requires a vertex with fit probability greater 
than 10^^, radius greater than 2 cm, and at most one hit per track inside of the vertex position. 
The x^ probability from the second method is shown in Fig. 
data and simulation. 



16c 



with good agreement between 



The nuclear interaction finder starts from the full list of tracks described in Section H] For 
each pair of tracks, the distance of closest approach is computed and if the two tracks are close 
enough they are considered linked together. A recursive finder produces blocks of tracks linked 
together from which a rough estimate of the displaced vertex position is computed. Finally, the 
tracks from a block are refitted together with a displaced vertex as a common constraint. V^ 
decays and photon conversions are removed from the resulting sample of displaced vertices. 
A tight selection is applied to the remaining vertices to remove fake tracks and pairs from the 
primary vertex. The resulting sample of significantly displaced vertices in the radial direction 
(r > 2.5 cm) is called the nuclear interactions sample. In the data, 80% of nuclear interactions 
are reconstructed with 2 tracks and 20% with 3 tracks. In the first case, a 30% combinatorial 
fake rate is expected from the simulation, while in the second case the fake rate is negligible. 



6.6 Study of b-tag Related Observables 
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Figure 16: Comparisons of data photon conversions (points) and real and fake photon conver- 
sion from simulation (filled histograms) showing: (a) distributions of the reconstructed pj of 
the converted photons from the first method, (b) the invariant mass of the e^e^ pairs from the 
first method, and (c) the distribution of the vertex x^ probability from the second method. The 
last bin of (b) is the overflow bin. 



The distribution of nuclear interaction positions provides a means of observing the material in 
the detector and validating the simulation of the material. The distribution of radial position r 



of the nuclear vertices, compared to the simulation, is shown in Fig. 17a The beam pipe at a 
radius of 3 cm, as well as the three barrel pixel layers at average radii of 4.3, 7.3, and 10.2 cm, are 
clearly seen. The radius is measured relative to the centre of the pixel detector. In the version 
of the simulation used here, this is also the centre of the beam pipe. In reality, the beam pipe 
centre is offset from the pixel detector centre resulting in a smeared distribution versus radius. 
Nevertheless, there is good agreement between the data and the simulation for the relative rate 
of nuclear interactions in the different barrel pixel structures and the beam pipe. This indicates 
a consistent description of the material distribution in this region. The material distribution 
in the endcap pixel detector is studied by selecting nuclear interactions with \z\ > 26cm and 
r < 19 cm. The longitudinal position |z| of the nuclear vertices, compared to the simulation, is 
show^n in Fig. 17b The pixel barrel flange (|z| < 30 cm) and the two pixel disks can be clearly 
distinguished. The tail up to 1 m is from pixel services. 

6.6 Study of Mag Related Observables 

The measurement of impact parameters and the reconstruction of secondary vertices have been 
tested with the limited event sample of December 2009. At higher collision energy these objects 
w^ill provide the main observables used in b-tagging algorithms. 

The 2009 data contain only a few w^ell-defined jets and mainly tracks at momenta below those 
typically used in b-tagging. To test the reconstruction on a sufficiently large sample, a few^ 
changes to the reconstruction chain have been applied with respect to what is described in 
Ref. [34|. As described in Ref. [35 1, jet reconstruction is performed using the anti-A:x jet cluster- 
ing algorithm II36II on objects obtained from the CMS particle flow reconstruction ||37ll38| . To 
recover low-momentum jets, the cone size is increased to 0.7 and the minimum pj is reduced 
to 3 GeV/c. The b-tagging algorithms are run on tracks associated with these jets. The track 
selection is also changed relative to Ref. I,34J : the minimum pj requirement is removed and 7 
rather than 8 hits are required. 

The impact parameter is computed with respect to the reconstructed primary vertex and the 
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Figure 17: Distributions of nuclear interaction vertices versus (a) radial position r for \z\ < 
26cm and (b) versus the magnitude of the longitudinal coordinate \z\ for \z\ > 26cm and 
r < 19 cm. The simulation histogram is normalized to the total number of nuclear interactions 
found in data in the full z range. 



distributions are compared between data and a minimum bias simulation reconstructed with 



the same algorithm settings. Figure 18a shows the three-dimensional impact parameter signif- 
icance distribution for all tracks in a jet. 

The secondary-vertex reconstruction using the tracks associated to jets has also been slightly 
modified compared to the algorithm described in Ref. [34] . The differences are a looser track se- 
lection, a relaxed vertex-to-jet direction compatibility, the use of track refitting in the secondary- 
vertex fit, and the use of the primary-vertex constraint rather than the beamspot. In addition, 
to suppress Kg candidates, the transverse secondary-vertex separation must be less than 2.5 cm 
and the secondary-vertex invariant mass more than 15 MeV/c^ from the nominal Kg mass. The 
significance of the distance between primary and secondary vertices is compared to what is ex- 
pected from a simulation of minimum bias events in Fig. |18b While many two- and three-track 
vertices are reconstructed, only one four-track vertex is found in the data. This event is shown 
in Fig. 19 



7 Conclusion 

The all-silicon CMS Tracker was designed for LHC physics. During its conception, emphasis 
was placed on providing fine granularity, large pseudorapidity coverage, and a number of 
redundant measurements to facilitate the pattern recognition. Meeting these design criteria 
has resulted in excellent tracking performance for the large-volume detector operating in a 3.8 T 
magnetic field in the early collision running. The proximity of the pixel detector to the LHC 
beam line permits precise reconstruction of primary and secondary vertices. The extended 
commissioning w^ith cosmic rays in 2008 and 2009 allowed most calibrations to be completed 
and provided a good initial alignment for most of the detector. This allowed fast and reliable 
operation of the Tracker in the first LHC collisions in December 2009. 

The performance of the Tracker has been studied using the collision data at centre-of-mass 
energies of 0.9 and 2.36 TeV. The tracking and vertexing resolutions are in agreement with the 



25 



Data 

Simulation (light jet) 
Simulation (charm jet) : 
Simulation (bottom jet) 




20 40 60 80 

signed 3D IP significance 



CD "-" 

O 

tr 

CD 

> 



CO 
T3 

o 
o 

CD 
CO 



10-^ 



1- 



Dala ' cms" 

Simulation (light jet) 
Simulation (charm jet) 
Simulation (bottom jet) 




5 10 15 20 25 30 

3D flight distance significance 



(a) (b) 

Figure 18: Distribution of (a) the significance of the three-dimensional impact parameter for all 
tracks in a jet and (b) the significance of the three-dimensional displacement of the secondary 
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Figure 19: Display of an event with a four-track secondary vertex. The vertex is separated from 
the primary vertex by 7(7 and the invariant mass of the four particles is 1.64GeV/c^, assuming 
they are all pions. 
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expected design performance in minimum bias events, for the level of the alignment achieved. 
Studies of the decays of Kg, A°, Er , and K*(892)^ test the capability to reconstruct displaced 
vertices and agree well with predictions from simulation. In particular, measurements of V^ 
mass, resolution, and lifetime provide strong tests of our understanding of the magnetic field. 
Tracker material, and detector performance. Our knowledge of the Tracker material is also 
evident in the agreement between data and simulation for photon conversions and nuclear 
interactions. Energy loss measurements in the Tracker, made possible by the calibration of the 
silicon strip analogue readout chain, provide good particle identification at low momentum as 
seen in the reconstructed (p and A decays. Finally, the alignment parameters determined in the 
cosmic ray studies are already adequate for the precise determination of impact parameters 
and the reconstruction of secondary vertices. This will ultimately be exploited for b-hadron 
physics and top-quark studies. 

In conclusion, the CMS Tracker has been commissioned to take data at the LHC. New collision 
data will allow more precise alignment and calibration, which will enable the study of the new 
vista of particle physics that is now opening up at the LHC. 
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